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Abstract 

The variational principle for Gibbs point processes with general finite range interaction is 
proved. Namely, the Gibbs point processes are identified as the minimizers of the free excess 
energy equals to the sum of the specific entropy and the mean energy. The interaction is 
very general and includes superstable pairwise potential, finite or infinite multibody potential, 
geometrical interaction, hardcore interaction. The only restrictive assumption involves the 
finite range property. 
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1 Introduction 


Gibbs point processes are popular models to describe the repartition of points or geometrical 
structures in space. They appeared first for modelling continuum interacting particles in statistical 
mechanics. Now they are widely used in as different domains as astronomy, biology, computer 
science, ecology, forestry, image analysis, materials science. The main reason is that they provide a 
clear interpretation of the interactio ns between the p o ints, such as attrac tion o r repulsion dependin; 
on their interdistance. We refer to IPrestoii ( IQTol . IChiu et ^ ( 2013ll . and Van Lieshoud ( 200' 
for classical text books on Gibbs point processes, including examples and applications. 

The Gibbs point processes are defined via their local unnormalized conditional densities of the 
form e~^ where H is an energy functional. They are the equilibrium states of the DLR equations 
(see definition 2). A variational principle, coming from the statistical physics, claims that the Gibbs 
measures are also the minimizers of the free excess energy equals to the entropy plus the mean 
energy. More precisely, for any stationary probability measure P on the space of configurations in 
the specific entropy X(P) with respect to the Poisson point process and the mean energy per 
unit volume H[P) are defined via thermodynamic limits (see ([S]) and ([7])). The principle claims 
that the Gibbs measures are exactly the probability measures P which minimize the functional 
P P{P) + H{P). It thus supports the common belief that the Gibbs measures provide a proper 
description of physical systems in thermodynamic equilibrium. They are many applications of the 
variational principle in physics and mathematics. In statistical mechanics, the phase transition 
phenomenon (non uniqueness of Gibbs measures) can be proved in studying the geometry of the 
set of Gibbs measures and in particu lar its extremal points. The variational principle is a key 
tool in this study (see iGeorgiil (|201l[l for a general presen tation). In pro bability theory, it is 
related to the large deviation principle for the empirical field Georgii ( 199411 . In spatial statistic, 
it is a crucial identifiabi l ity as sumption for the consistency of the maximum likelihood estimator 
Dereudre and Lavancied ( 2015ll . This last recent paper highlights the importance of the variational 
principle for models coming from the spatial statistic. It was our initial motivation for the present 
paper. 
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For the lattice Gibbs m odels, the variational principle is well established an d a general proof can 
be fou nd in IPrestoM ( 1976 ). Section 7. The first proof was for the Ising mo del in Landlord and Ruellel 
( 1969ll . In the setting of Gibbs point processes, they are less results. In iGeorri il ()l994bll3 l the au¬ 
thor proves the variational principle for the pairwise potential energy H{uj) = — y\) 

where the sum is over all couples {x,y} in the configuration ui. The potential cj) is assumed to be 
non-integrably divergent at the origin (i.e. ds = oo) producing a strong repulsion when 

the particles are closed to each other. A typical examples is the Lennard Jones pairwise potential 
4>(r) = —br~^. In another wo rk the variational principle is proved for the Delaunay-tile inter¬ 

action Dereudre and Georsfiil ( 2009| i. The energy function has the following form H{uj) = 
where the sum is over all triangles T of the Delaunay triangulation based on w. It is a contin¬ 
uum spatial version of nearest neighbours interaction models. As fas as we know both papers are 
the only ones proving the variational principle for Gibbs point processes models. Unfortunately 
many interesting energy functions are not covered by these results, as for example any pairwise 
potential energy with bounded potential (j). In particular, the well-studied Strauss model in spatial 
statistics with the pairwise potential (f>(r) = I1[o,k] (r) is uncovered. In s tochastic geometry, the 
Area- inter action or the Querma ss-interaction are not covered as well (see iBaddelev and Lieshout 
( 1995 ) and Kendall et al. ( 19991) ). 

In this paper we prove the variational principle for Gibbs point processes with general finite 
range interaction. The other assumptions are standard and satisfied by all models we met in 
statistical mechanics and spatial statistic (see Section 3). In particular, our setting includes super¬ 
stable pairwise potential, finite or infinite multibody potential, geometrical interaction, hardcore 
interaction. The proof is based on fine controls of the relative entropy of Pa with respect to the 
Gibbs measure on A, where P is any statio nary field on and A an obser vable window tendin g 
to This strategy was already present in Preston ( 19761) . Georgii ( 1994a ) and Georgiil ( 1994bl) . 

The paper is organized as follows. In section 2, we introduce the notations and the Gibbs 
models. The variational principle and the main theorem are presented in Section 3. Two standard 
examples are given in Section 4; the superstable pairwise potential with compact support and the 
Quermass interaction. Section 4 is devoted to the proof of our main theorem. 


2 The Gibbs models 

2.1 State spaces and reference measnres 

Our setting is the Euclidean space of arbitrary dimension d > 1 equipped with its Borel tr- 
field. An element of is denoted by x and the Lebesgue measure on is denoted by X‘^. A 
configuration is a subset ui of which is locally finite, meaning that w O A has finite cardinality 
Na{uj) = #(w n A) for every bounded Borel set A. The space D of all configurations is equipped 
with the (T-algebra P generated by the counting variables Na- The space of finite configurations 
is denoted by D/. 

The symbol A will always refer to a bounded Borel set in R'^. It will often be convenient to 
write ujA in place of w D A. We abbreviate lu U {a;} to w U a; and abbreviate w\{a;} to U!\x for every 
w and every a; in w. 

As usual, we take the reference measure on (D, to be the distribution tt of the Poisson point 
process with intensity measure on R*^. Recall that tt is the unique probability measure on (D, P) 
such that the following hold for all subsets A: (i) Na is Poisson distributed with parameter A‘^(A), 
and (ii) conditional on Na = n, the n points in A are independent with uniform distribution on A. 
The Poisson point process restricted to A will be denoted tta. 

Translation by a vector u S R'^ is denoted by t^, either acting on R'^ or on D. A probability P 
on D is said stationary if P = P o for any u in R'^. In this paper we consider only stationary 
probability measures P with finite finite intensity measure (i.e. Ep(A^[q ijd) < -|-oo). We denote 
by P the space of such probability measures. 
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2.2 Gibbs point processes models 

We consider a measurable function H from Q.f to K. U {+ 00 } which is called energy function. We 
assume that H is stationary; for any uj € ^If and any u G = H{tu{uj). We assume 

also that the energy function H is hereditary which means that for any x in and uj in ilf, 
H{ujU{x}) = +00 as soon as H{uj) = + 00 . The energy H is said non-degenerate if iJ({0}) 7 ^ 00 
and H{%) = 0. We assume also that H is stable which means that there exists d > 0 such that 
for any finite configuration uj G flf 


H{uj) > -AN{uj). (1) 

All these assumptions are standard and non restrictive. The main restriction in the present 
paper is the finite range assumption which means that there exists R > 0 such that for any 
configuration uj, any bounded set A the quantity 


Ha{uj)-.= H{uja')-H{ujj,,\^) (2) 

(with the convention 00 — 00 = 0) does not depend on the choice of A' as soon as A©i3(0, R) C A'. 
Ha{uj) represents the energy of uja inside A given the configuration outside A. 

The Gibbs measures P associated to H are defined through their local conditional specification, 
as described below. We denote by floo the set of configurations uj G ^ such that for any A, 
H(uja) < +00. So for every A and every configuration w G floo, the local conditional density /a of 
P with respect to tta is defined by 


/a(w) 


■^a(wa‘=) 


(3) 


where .^a(wa<=) is the normalization constant given by 


2'a(wa=) = J ^^^‘^^^‘^^‘=^7rA(da;A). 


Let us note that 0 < Za(wa<=) < +00 since H is stable and non-degenerate. 

We are now in position to define the Gibbs measures associated to H fSee lPreston ( 1976ll for 
instance). 


Definition 1 A probability measure P on is a Gibbs measure for the energy function H if 
P(Ooo) = 1 and if for every bounded borel set A, for any measurable and bounded function g from 
fl to R, 

J g{uj)P{duj) = J j ^(wa Ua;A<=)/A(wA U WAO’^A(da;A)P(da;). (4) 

Equivalently, for P-almost every uj the conditional law of P given uja‘= is absolutely continuous with 
respect to tta with the density /a defined in ©. 


The equations (|3]) are called the Dobrushin-Lanford-Ruelle (DLR) equations. The existence 
of such Gibbs measu res, in the present setting of finite range stable interactions, is done in 
Dereudre et al.l ( 2012ll . Corollary 3.4 and Remark 3.1. Note that the uniqueness of such P does 


not necessarily hold, leading to the phase transition phenomenon. We denote by Qh the set of all 
Gibbs measures for the energy H. 


3 Variational Principle 

The variational principle in statistical mechanics claims that the Gibbs measures are the minimizers 
of the free excess energy dehned by the sum of the the mean energy and the specific entropy. 
Moreover the minimum is equal to minus the pressure. Let us first define precisely all these 
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macroscopic quantities. For the sake of simplicity we consider the macroscopic limit along the 
sequence of sets An = [—n,n]‘^, n > 1. Limits A in the Van-Hove sens could have been 

considered as well. 

Let P be a stationary probability measure in V. The specific entropy of P is defined as the 
limit 


I(P) = lirn -^X(PA„,7rA„), 

n-H-oo |A„| 

where for any probability measures /i and v 


(5) 




/ ln(/)d/r ii fj, ly with density / 
- 1-00 otherwise 


Note that the limit in JS]) always exists; see Georgii ( 201lll for general results on specific entropy. 
Let us now introduce the Pressure. It is defined as the following limit. 


PH 


:= lim 

n—¥-\-oo 


\An 


■ln(Z„), 


( 6 ) 


where = Za„ (0 ) is the partition function with empty boundary condition. 

In the following lemma we show that pn always exists in the setting of the present paper. 


Lemma 1 Assuming that the energy function H is finite range, stable and non-degenerate, then 
the pressure pn defined in (O exists and belongs to [—1, (e^ — 1)]. 

Proof. 

For any set A we denote by A® the set 

A® = {a; G A, B{x, Pq) C A}, 

where Pq is an integer larger than the range of the interaction P. So for n > Pq, A® = A„_/Jq. 
For any Pq < m < n, we consider the Euclidean division n = km + I with 0 < I < m, k > 0. 
Let {Aln)i<i<k‘‘ a family of k‘^ disjoint cubes inside A„ where each cube is a translation of Am- 
From the definition of the partition function 




An\UiA: 


= 0) y ( 




'UiA. 


:,e{duj) 




II^« 


— ® ^m-Ro' 

So since \An\/k'^ goes to |Am| when n goes to infinity, 

liminfln(Z„) > -^ {Zm-Ro - 2dRo{2m)‘^~^). 
"->■00 |A„| |Am| 

This inequality holds for each m > Rq. So, letting m tends to infinity 

lim inf - ln(Z„) > lim sup -j-—= limsup -j-— tZ„ 

n^oo |A„| m->-oo |A m\ m—¥oo IAttt,! 

which proves that the limit exists in M U {±oo}. 

Thanks to the stability and the non degeneracy of H we get that 
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which implies that pn S [—1, (e^ — 1)]. 

■ 

The last macroscopic quantity involves the mean energy of a stationary probability measure 
P. It is also defined by a limit but, in opposition to the other macroscopic quantities, we have to 
assume that it exists. The proof of such existence is based on stationary arguments and a nice 
representation of the energy contribution per unit volume. Examples are given in Section |4l So 
for any stationary probability measure P we assume that the following limit exists 

H{P) := lim f H{uJAjPiduj). 

n^oa |A„| J 

and we call the limit mean energy of P. 

We need to introduce a last technical assumption on the boundary effects of H 
that for any P in Qh 

lim [ dHA^{u})P{(hj) = 0, (8) 

rt^oo |A„| J 

where = Ha„{uj) — H{uja„)- This assumption is satisfied by all the examples we meet. 

Theorem 1 We assume that H is stationary, hereditary, non-degenerate, stable and finite range. 
Moreover we assume that the mean energy exists for any stationary probability measure P (i.e. 
the limit ([7]) exists) and that the boundary effects assumption ([8]) holds. Then for any stationary 
probability measure P € P 


(7) 

. We assume 


I{P) + H{P) > -PH, 

with equality if and only if P is a Gibbs measure (i.e. P G Gh)- 


(9) 


4 Examples 

In this section we present two examples of energy functions included in the setting of Theorem 
[TJ The first example is the standard superstable pairwise potential energy. The second example 
involves the Quermass interaction which is an energy function for morphological patterns built by 
unions of random convex sets. It can be also viewed as a infinite body potential interaction. The 
main restriction in Theorem [T] is the finite range property and so any standard examples, having 
this property, could have been considered as well. 


4.1 Pairwise potential 

In this section the energy function has the following expression: for any finite configuration uj G Gif 


H{uj) = zN{uj) + ^ (j){x-y), 


( 10 ) 


{x,v}Cuj 


where is a symmetric function from to M U {+oo} with compact support. The parameter 
z > 0 is called activity and allow to change the intensity of the reference Poisson point process. 
The potential (p is said stable if the associated energy H in (I10|) is stable. In the following we 
need that the potential (p is superstable which means that <p is the sum of stable potential and a 
positive potential which is non negative around the origin. See Ruellel ( 1969ll for examples of stable 
and superstable pairwise potentials. In this setting the variational principle holds as a corollary of 
Theorem [TJ 
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Corollary 1 Let H be a energy function coming from a superstable pairwise potential (j) given by 
(uni). Then for any stationary probability measure P € V 


I{P) + H{P)>-ph. (11) 

with equality if and only if P is a Gibbs measure (i.e. P € Gh)- The expression of H{P) is 
given in m- 

Proof. 

Let us check the assumptions of Theorem [T] It is obvious that H is stationar y, hereditary no n 
degenerate, stable and hnite range. The existence of the mean energy is proved in ICeorgiil (jl994br i. 
Theorem 1. It is given by 


H{P) = 


i f Eo^xeuj </>(x)P^(du;) if < oo 

+00 otherwise 


( 12 ) 


where is the Palm measure of P. Recall that can b e viewed as the natural version of the 
conditional probability P(.|0 £ w) Isee iMatthes et al~ ( 1978h for more details). So, it remains to 
prove the boundary assumption ([5]). Let P a Gibbs measure in Qp. A simple computation gives 
that for any w £ 


where A® = An+Rg and A® = An -with Ro an integer larger than the range of the interaction 
R. Using the GNZ equation (see Nguven and Zessin ( 1979lP . the stationarity of P we obtain 


\EpidHAj\ < 


Y Y -2/)l^(c^w) 


H^-V) 

An\A,l 


g - yj 1 ^( 3 ; _ y)\dxP{ckjj) 


|A®\A„|e ^ J e ^ \<p{y)\P{duj). 


yei^B(o,Ro) 


Since (f is stable we deduce that (j) > —A—2z. So denoting by C := sup 2 g[_^_ 2 z;_|_oo) |c|e '^ < 00 
we find that 


\Epmf^E\ < |A®\A„|Ce-^y'fVB(o.K„)(cc)e(-^+"^)'^-®.«o)Mp(da;). (13) 

Using the estimates in Ruelle ( 1970ll corollary 2.9, the integral in the right term of (fT^ is finite. 
The boundary assumption ([5]) follows. 


4.2 Quermass interaction 

The Quermass process is a morphological interacting model introduced in Kendall et al.l (19^ _ 

which i s a generalization of the well-kn own Widom-Rowlinson process or Area Process Isee lWidom and Rowlinson 
(197^, Baddelev and Lieshoud ( 1995ll '). Since the existence of the Quermass process is only proved 
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in the case d < 2 we restrict the following to the non trivial case d = 2. For any finite configura¬ 
tion w, L{uj) denotes the set \Jx^cjB{x,r) and the energy is defined as a linear combination of the 
Minkowski functionals; 


H{uj) — 9iA(^L{uj)j + 02£-(^L{uj)j + 03 x(^L{uj)^ , 


(14) 


where r > 0, Oi GR,i = 1...3 are parameters and A, y are respectively the area, the perimeter 
and the Euler-Poincare characteristic functionals. Recall that x(L(w)) is equal to Ncc{L{uj)) — 
Nh{L{u))) where Ncc{L{uj)) deno tes the number of connected components in L{uj) and Nh{L{uj)) 
the number of holes. We refer to IChiu et al.l (j2013|) for more details about Minkowski functionals. 


Corollary 2 Let H be the Quermass interaetion given in m- Then for any stationary probability 
measure P € V 


I{P)+H{P)>-ph, (15) 

with equality if and only if P is a Gibbs measure (Le. P € Qh)- The expression of H{P) is 
given in (HU). 

Proof. 

As in the previous section, we check the assumptions of Theorem [TJ It is obvious that H 
is stationary, hereditary and non degenerate. In dimension d = 2 the functional x satisfied the 
following bound 


|x(L(cc))| <3iV(a;), 


(16) 


see 


Kendall et al.l 1 1999ll . The stability of H follows easily. The finite range assumption is a 


consequence of the additivity of Minkowski functionals. Note that the range of the interaction is 
R = 2r. In the following we denote by C the cube [0,1]^, by dC the boundary of C and by C the 
double edges {0} x [0,1] U [0,1] x {0}. For any k & I? we consider also the translations Ck = Tk{C), 
dCk = Tk (dC) and Ck = Tk (C). Thanks to the additivity of Minkowski functionals we obtain that 
for any finite configuration uj and any n > 1 


h{cvaJ = Yi (eiA{L{io)nCk)+e2 c{L{cv)nCk)-c{L{oj)ndCk) 

kG{—n,n — l}^ 


+03 


x{l{lo) n Ck) - Ncc{l{lo) n Ck) ] -i- i?„(wA„), 


(17) 


which gives the energy contribution of each cube Ck in H{uja„)- Thanks to (fTOl) and obvious 
bounds for A and £, the boundary term i?„(a;A„) satisfies for some constant c > 0 

|i?n(wA„)| < cN (u;a„\a„_bo) ■ 

For any stationary probability measure P £ P we deduce easily from (HID and ([T^ the existence 
of the mean energy © with 


H{P) = / 0iA{L{u;)nC)+02[c{L{uj)nC)-C{L{u;)ndC)\+03[x{Li^)<^C)-Ncc{L{io)nC))P{du;). 

(19) 

Thanks to (HZl and m, the boundary assumption ([5]) is satisfied as well. 
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5 Proof of Theorem [T] 


Let us start by proving the inequality ([5]) for any stationary probability measure P S P. For u > 1 
we define the Gibbs measure on A„ with free boundary condition by 


QniduJAj = -^e 


So 


APa„,Qu) = 


In ( ^^(a;A„) ) dPA„(wA„) 


In 


dQu 

dPKr. 


(‘^A„)^^^^(wa„) ) (iPA„(wA„) 


d7rA„ " dQr, 

= 2:(PA„,7rA„) + J H{uJAjdP{u;) +\n{Zn), 


which implies that 


71“ 

n^oo |A^ 


{Pa„ , Qn) = AP) + H{P) + PH. 


( 20 ) 


( 21 ) 


Since Ia„ (Pa„ j Qn) is positive the inequality dH) follows. Let us now prove that for any P & Gh 
the equality holds in © . Let us show that the limit in (I^Tl) is negative. Recall that Rq is an integer 
larger than the range of the interaction R and that A® stands for the set An+Ro- We denote by 
'A'An G ^A®\A„ point process on A® with independent configurations on A„ and 

A®\A„ with distributions tta^ and Pa®\a„ respectively. Then 


lim I(PA„,Qn) = lim 
n->oo |A„| “ 


1 


= lim 


ra->oo |A„ 
1 


, Qn+Ro ) 

dP^Qf 


\A„ 


In 


d-KAn G^a®\a, 


dnAn^PAf\Anr , 
- Aa^ )--(‘^A®) 


A® 


d-TT. 


dQ 


n+Ro 


(Wa®) dA’A®(wA®) 


(^(-^A®\A„>^A®\A„) + MZn+Ro) + J - ^A®(^A®) 

— ln(^A„(w^®))dP^®(w^®)y (22) 


where the densities w hich appear abo ve are given by (|3]) and (1^ . By subadditivity of the entropy 
( Proposition 15.10 in Georgii ( 201lll l. 


0 < P{Pa® \A„ ’ ’’’a? \A„ ) — ^(Pa^ \A„ ’ ’’’a? \A„ ) dl(PA„ , TTAn ), 

which implies that 


|A„y^'^^®\A"’’"A®\Aj - 0- 

Moreover, thanks to the existence of the mean energy and the boundary assumption assumption 
dH), the term lim„_>oo f(H(ujj^e) — (uj)dP(uj) vanishes as well. Therefore the limit in 

dSD is negative provided we show 





















liminf f In dP{uj) > 0. (23) 

n->-oo \An\ J V -^n+Ro J 

From the definition of Z\^{u!), 

Za„{uj)> [ Hr , 1 

and therefore 


lim inf -— 

n-s-oo |A„ 


In 


Za^ (^) 




TI-\-Rq 


dP{uj) > liminf —- (ln(Z„__Rj - \n{Zn+Ro) - |A„\A„__rJ) = 0 

n^co |A„| 


which proves (1^^ . The proof of Theorem [T] is complete if we show that any stationary probability 
measure P sol ving the equality in ([5]) is a Gibb s measu re. We follow essentially the scheme of 
IPrestonI ( 1976ll (In the variant used in iGeorgii (1994a), Section 7). So let P be a stationary 
probability measure such that X{P) + H{P) + pn = 0. Let us show that for any bounded local 
function g and any bounded set A, f g(u})P(duj) = f gA{uj)P{uj) where the function gA is defined 

by 


ffA(w) = y g(wA Ua;A<=)/A(wA Ua;A<=)7^A(da;A). 


Without loss of generality we assume in the following that \g\ is bounded by one. Thanks to 
the equality (1^ . for n large enough I(Pa„, Qn) is finite and therefore Pa„ admits a density with 
respect to Qn which we denote by /„. Let A' be a bounded set such that A® C A' and such that 
g is Pa'- measurable. For n large enough such that A' C A„, the probability measure Pa' admits a 
density with respect to Qn restri cted to A ' which we denote by fn,A'- Since X(Pa„, Qrt)/|A„| —>• 0, 
using the standard Lemma 7.5 in lGeorgiil ( 1994a ). for any <5 > 0 there exists n large enough and a 
set A' with A® c A' c An such that 


J \fn,A' /ri,A'\AMQ™ ^ 


We obtain that 


J g(w)-5A(w)P(dw) = y /„,a'(‘*^)5(w) -/„,A'\A(w)gA(w)(5n(dw)■ 

From the definition of Qn and since A® C A„ we have 

y fn,A'\-^(^)9A{^)Qn{duj) = J /„,A'\A(w)g(a;)Q„(dw) 

and we deduce that | / g{uj) — gA{oj)P{duj)\ < 6. Letting S tends to zero we get the DLR equation 
on A. The proof of Theorem [T] is complete. 
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